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Although galactic dark matter halos are basically Newtonian structures, the study of their in- 
terplay with large scale cosmic evolution and with relativistic effects, such as gravitational lenses, 
quintessence sources or gravitational waves, makes it necessary to obtain adequate relativistic de- 
scriptions for these self-gravitating systems. With this purpose in mind, we construct a post- 
Newtonian fluid framework for the "Navarro-Frenk- White" (NFW) models of galactic halos that 
follow from N-body numerical simulations. Since these simulations are unable to resolve regions very 
near the halo center, the extrapolation of the fitting formula leads to a spherically averaged "univer- 
sal" density profile that diverges at the origin. We remove this inconvenient feature by replacing a 
small central region of the NFW halo with an interior Schwarzschild solution with constant density, 
continuously matched to the remaining NFW spacetime. A model of a single halo, as an isolated ob- 
ject with finite mass, follows by smoothly matching the NFW spacetime to a Schwarzschild vacuum 
exterior along the virial radius, the physical "cut-off" customarily imposed, as the mass associ- 
ated with NFW profiles diverges asymptotically. Numerical simulations assume weakly interacting 
collisionless particles, hence we suggest that NFW halos approximately satisfy an "ideal gas" type 
of equation of state, where mass-density is the dominant rest-mass contribution to matter-energy, 
with the internal energy contribution associated with an anisotropic kinetic pressure. We show that, 
outside the central core, this pressure and the mass density roughly satisfy a polytropic relation. 
Since stellar polytropes are the equilibrium configurations in Tsallis' non-extensive formalism of 
Statistical Mechanics, we argue that NFW halos might provide a rough empirical estimate of the 
free parameter q of Tsallis' formalism. 

PACS numbers: 04.20.-q, 02.40.-k 



I. INTRODUCTION 

A large amount of compelling evidence based on di- 
rect and indirect observations: rotation velocity profiles, 
microlensing and tidal effects affecting satellite galaxies 
and galaxies within galaxy clusters, reveals that most 
of the matter content of galactic systems is made up of 
dark matter (DM). Since the physical nature of DM so 
far remains uncertain, this issue has become one of the 
most interesting open problems in astrophysics and cos- 
mology P, 0, La, 0, Ei- Among a wide variety of pro- 
posed explanations we have: thermal sources, meaning 
a colissionless gas of weakly interacting massive particles 
(WIMP's), which can be very massive (to ~ 100 — 200 
GeV) and supersymmetric Q ( "cold dark matter" CDM) 
or self-interacting less massive (m ~ KeV) particles 0)0 
("warm" DM). 

The DM contribution dispersed in galactic halos is 
about 90-95 % of the matter content of galactic systems, 
while visible baryonic matter (stars and gas) is clustered 
in galactic disks. It is then a good approximation to 
consider the gravitational field of a galaxy as that of its 
DM halo (for whatever assumptions we might make on 
its physical nature), while visible matter can be thought 
of as "test particles" in this field 0, . 

Assuming the CDM paradigm, we can distinguish two 
types of halo models: idealized models obtained from a 



Kinetic Theory approach, whether based on specific the- 
oretical considerations or on convenient ansatzes that fix 
a distribution function satisfying Vlassov's equation , 
or models based on "universal" mass density profiles ob- 
tained empirically from the outcome of N-body numeri- 
cal simulations [ll|, [l^, LU| . In this paper we will study 
the equilibrium configurations that emerge from the lat- 
ter approach, based on the well known numerical sim- 
ulations of Navarro, Frenk and White (NFW) [U lilf . 
It is important to mention that these simulations yield 
virializcd equilibrium structures that reasonably fit CDM 
structures at a cosmological scale (> 100 Mpc), though 
some of their predictions in smaller scales ( "cuspy" den- 
sity profiles and excess substructure) seem to be at odds 
with observations [TEl Il6j . especially those based on 
galaxies with low surface brightness (LSB), which are 
supposed to be overwhelmingly dominated by DM and 
so well suited to examine the predictions of various DM 
models [HEl. 

We consider in this paper that DM halos are spheri- 
cally symmetric equilibrium configurations, a reasonable 
approximation since their global rotation is not dynami- 
cally significant Galactic halos in virialized equilib- 
rium are also Newtonian systems characterized by typical 
velocities, ranging from a few km/sec for dwarf galaxies 
up to about 1000 — 3000 km/sec for rich clusters. So, why 
bother with a general relativistic treatment? First, there 
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is a purely theoretical interest in incorporating these im- 
portant self-gravitating systems into General Relativity, 
the best available gravitational theory. In fact, impor- 
tant experimental tests of General Relativity are cur- 
rently and customarily carried on (within a weak field 
post-Newtonian approach) in Solar System bound New- 
tonian systems. Secondly, a post-Newtonian description 
of galactic halos can be, not only useful and interesting, 
but essential for studying their interaction with physical 
effects that lack a Newtonian equivalent, such as gravi- 
tational lenses or gravitational waves. In fact, the post- 
Newtonian halo models that we present in this paper 
can be readily used in lensing studies, or can provide 
the unperturbed zero order configuration in the exami- 
nation of the perturbative effect of gravitational waves 
on galactic halos. Also, a post-Newtonian description is 
necessary in the study of the interplay between galactic 
structures and large scale (> 100 Mpc) cosmological evo- 
lution dominated by a repulsive "dark energy" , modeled 
by sources like quintessence and/or a cosmological con- 
stant, whose Newtonian description might be inadequate. 
Finaly, since galactic halos are customarily examined as 
Newtonian structures, we feel it is important to show 
the readership of General Relativity journals how to con- 
struct spacetimes, in a post-Newtonian approximation, 
that are suitable for the description and study of these 
important self-gravitating systems. 

Given the NFW mass density profile, we show in sec- 
tion II how all dynamical variables of the Newtonian 
NFW halo can be derived. In section III we construct 
a post-Newtonian fluid relativistic generalization of a 
NFW halo, under the assumption that the gas of col- 
lisionless WIMPs should satisfy an "ideal gas" type of 
equation of state pol l2ll . For isotropic velocity dis- 
tributions, this assumption allows us to determine the 
internal energy density by means of the hydrostatic equa- 
tions themselves. For the case with anisotropic velocities 
we follow the same procedure with regards to the "radial" 
component of the stress tensor, determining the "tangen- 
tial" stress by a suitable empirical ansatz (section VI-B). 

The fact that the NFW spherically averaged mass den- 
sity profile diverges at the halo center follows from in- 
terpolating a fitting formula associated with numerical 
simulations that have a finite resolution limit at the halo 
center j23| ■ Although this behavior of the density profile 
does not imply that simulations predict an infinite cen- 
tral density, it is none-the-less an undesirable feature, 
which we ammend in section IV by replacing a small 
region around the center of the NFW halo with a spher- 
ical section of a spacetime with constant matter-energy 
density, i.e. a section of the "interior" Schwarzschild so- 
lution, that is continuously matched to the remaining of 
the NFW post-Newtonian spacetime. 

Galactic halos are hierarchical structures: small halos 
lie within galaxy clusters, which might be part of su- 
perclusters, etc, the asymptotic field of a typical NFW 
halo should somehow merge with a mean field of a larger 
substructure, or with a mean cosmological field. How- 



ever, the dynamical input from a suitable cosmological 
background is well imprinted in the theoretical design of 

NFW simulations [HE El, while the effect of a cos " 
mological constant on the equilibrium of virialized halo 
structures is know to be negligible (see [24| and refer- 
ences quoted therein). Hence, at galactic scales the em- 
piric NFW profiles are assumed to be valid only up to 
the "virial radius" , a physical "cut-off" scale associated 
with a virialization process 0, 0, El El E3 > ignoring al- 
together their transition to background fields associated 
with larger structures or to a cosmological background. 
While a post-Newtonian approach also allows one to im- 
pose this virial cut-off scale and to ignore its asymptotic 
behavior, we show in section IV how well behaved asymp- 
totically flat NFW configurations can be constructed. 
Also, since any localized self-gravitating system (even if 
belonging to large substructures) can be approximately 
described as an isolated system, we also examine an al- 
ternative cut-off by matching generic NFW halos to a 
Schwarzschild vacuum exterior at the virial radius. 

In section V we provide the equilibrium equations given 
in terms of suitable dimensionless variables for the post- 
Newtonian NFW halos, using the matching with the "in- 
terior" and "exterior" Schwarzschild solutions defined in 
section IV. Analytic solutions of these equations are ob- 
tained in section VI, for isotropic velocities (subsection 
A) and for a well defined case with anisotropic velocities 
(subsection B). We show in both cases that (outside the 
central region) the radial pressure and mass density sat- 
isfy approximately a polytropic relation characteristic of 
stellar polytropes [l(j . Even if NFW halos exhibit (in 
general) deviations from an isotropic velocity distribu- 
tion, while velocities in polytropes are strictly isotropic, 
we argue in section VII that the resemblance of outer 
regions of NFW halos to stellar polytropes might be sig- 
nificant, since virialized self-gravitating systems exhibit 
non-extensive forms of energy and entropy, and stelar 
polytropes are the equilibrium states in the application 
to astrophysical systems of the non-extensive Statistical 
Mechanics formalism developed by Tsallis |S EE E3. 
(see [2^ for a critical approach to this formalism). 



II. THE NFW DARK MATTER HALOS. 



The well known N-body numerical simulations by 
Navarro, Frenk and White (NFW) yield the following 
"universal" expression for the density profile of virialized 
galactic halo structures [IJ [II El [3 



x(l + xf ' 
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while the concentration paramet er en can be expressed in 
terms of the virial mass M vi[ by [30( 



c = 62.1 x 



(M wit h 



-0.06 



(5) 



where —1/3 < e < 1/2. The virial radius r vir is given in 
terms of M vir by the condition that average halo density 
equals A times the cosmological density po 



A Po = 



3 M vir 



(6) 



where A is a model-dependent numerical factor (for a 
ACDM model with total Q = 1 we have A ~ 100 31]). 
Hence all quantities depend on a single free parameter 
M vir with a dispersion range given by e for different halo 
concentrations . 

Using this profile, the mass function and gravitational 
potential follow from the Newtonian equations of hydro- 
static equilibrium 



M' = inpr 2 , 
*/ GM 



(7) 
(8) 



where a prime denotes derivative with respect to r. 
Hence, the NFW mass function follows from integrating 
@ for p given by (JIJ 



M N 



4 7r r 3 <5 Po 



ln(l + x) 



1 



(9) 



so that M NFW (0) = 0, while M NFW evaluated at r = r vir 
(or x = Co) yields M vir as defined in (JHJ for 8q given 
by Circular rotation velocity and the gravitational 
potential follow from © 



V' 1 



= -V? 



ln(l + x) 

x 

ln(l + x) 



1 



where the characteristic velocity Vq is 



y 2 = inGr 2 s S oPo = 



_ 3S Q GM vir 



(10) 
(11) 

(12) 



and the integration constant was chosen so that # NFW — > 
as x — > oo. Notice that, even if p NFW diverges, all 



other quantities (but not their gradients) are regular as 
as x — > 0. 

Since numerical simulations usually yield anisotropic 
velocity distributions, we have in general an anisotropic 
stress tensor so that "radial" and "tangential" pressures, 
P = P r and P± are involved in the Navier-Stokes equa- 
tion 
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P'= - P 9'-—P, 



where 



r = 



P-Pa 



P 



(13) 



(14) 



is the anisotropy factor. Given p NFW and M NFW , the ra- 
dial and tangential pressures follow from integrating i|13|) 
for a given choice of T. For the NFW forms JIJ and J^J, 
there are analytic solutions of (|13fl for T = (isotropic 
case) and for various empiric forms of T |l4j . 



III. RELATIVISTIC GENERALIZATION 

Under the assumptions that we outlined in the Intro- 
duction, the spacetime metric for an NFW dark matter 
galactic halo should be a particular case of the spherically 
symmetric static line element 



ds 2 



< xi. ( ) c 2 dt 2 



1 - 



2GM 



dr 2 



r 

+r 2 {d9 2 + sin 2 6> d</> 2 ),(15) 



so that M(r) has units of mass. The functions $(r) and 
M (r) are suitable relativistic generalization of the NFW 
functions given by © and iJTTJ. We will assume a fluid 
energy-momentum tensor of the most general form for 
the metric (|15fl 



T 



ab 



nu a u b + P h a0 + n a0 , 



ab I Tjab 



(16) 



where p and p are the matter-energy density and 
isotropic pressure along a 4-velocity field u a = 
exp(-$/c 2 ) S a t , while h ab = g ab + u a u b and n afc is the 
anisotropic and traceless (n a a = 0) stress tensor, which 
for the metric (|15|) takes the form 



n a 6 = diag [o, -2n, n, n] 



(17) 



so that p and n = H(r) relate to the radial and tangential 
pressures, P and P±, by 



P±-P = 3n, 2P X +P = 3p. 



(18) 



The field equations and momentum balance (T ab ; j = 0) 
associated with (|15H - <|18H are 



M' = 4irpr 2 /c 2 , 

G [M + 4 7r P r 3 / c 2 ) 



P' 



r[r-2GM/c 2 ] ' 

N $' 2T 
(p + P)^ P, 



(19) 
(20) 

(21) 
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where T is given by l|14|l . These equations are the 
relativistic generalization of the Newtonian equilibrium 
equations JJJl, JHJ and l|13(l . In the Newtonian case all 
these equations are decoupled, so that once p is known 
and r is prescribed, all other quantities follow by sim- 
ple integration of quadratures. In the relativistic case 
we have, in general, three equations for five unknowns 
{p. P, T, M, $). Thus, we must provide a relation be- 
tween p and P, together with a suitable assumption that 
determines or prescribes the form of T. 

Since the WIMPs in the collisionless gas making up 
galactic halos are interacting very weakly, it is reasonable 
to consider such a gas as approximately an "ideal gas" 
whose total matter-energy density, p,, is the sum of a 
dominant contribution from rest-mass density, p c 2 , and 
an internal energy term that is proportional to the pres- 
sure P and to the velocity dispersion a 2 = (v 2 ) ~ (v±)- 
Hence, we shall assume that the matter source of NFW 
halos complies with the equation of state of a non- 
relativistic (but non- Newtonian) ideal gas |1 El I21I E^] 



pc 



3 a 2 
2 c 2 



P 



pa 



(22) 



where we can identify p = /3 NFW (or with any mass density 
formula used in halo models) and the velocity dispersion 
is related to a kinetic temperature T by 01 



P 
P 



(23) 



where k B is Boltzmann's constant. At this point, we 
believe it is convenient to mention the following two ide- 
alized models of self-gravitating systems as useful theo- 
retical references liol: 



Isothermal Sphere: 

? kr. To 
cr 2 = -2— - = const, 
m 

Stellar Polytropes: 

2 j-x I In 
(7 = K p 1 => 



P = Kp, 



P = Kp 1 



-l/n 



(24) 



where To, K and n (polytropic index) are constants. The 
isothermal sphere corresponds to a Maxwell-Boltzmann 
velocity distribution, the equilibrium state associated 
with the extensive Boltzmann-Gibbs entropy [Toll2lLl22]| . 
The stellar polytropes are also solutions of the Vlassov 
equation, but are associated with the equilibrium state in 
the non-extensive entropy functional proposed by Tsal- 
lis IMEE IH HH. Notice that the isothermal sphere 
follows from the stellar polytropes in the limit n — > 00 
(the extensivity limit in Tsallis' formalism). 

For Newtonian characteristic velocities in galactic ha- 
los, we have tr 2 /c 2 <C 1 and p » pc 2 and so P ~ P± <C 
pc 2 , so that IL'l'l) provides a plausible equation of state 
for a relativistic generalization of galactic halos. It is ev- 
ident that in the Newtonian limit a 2 j c 2 — > we recover 



the Newtonian equilibrium equations (0, JHJ and 11311 . 
What needs to be done now is to insert the equation of 
state (1221) into the field equations 1191) (J2TJ. It turns out 
to be easier to work with P instead of a or T, using l)22|) 

as 



pc 



(25) 



so that a and/or T can be obtained afterwards from P 
through J32Jl and J33J). Combining and (JSIJ into a 
single equation and using 1)25(1 we obtain the set 



M' 
P' 



= 47T 



3 P 
2 

G[p + \Pjc 2 \ [M + 4-7T (P/c 2 ) r 3 ] 2T 



(2G/c 2 ) M] 



(26) 



P. 

(27) 



which becomes determined once we identify p = p NFW 
and specify Y — Y(r). We can solve these equations in 
a post-Newtonian approximation by keeping only terms 
up to order a 2 /c 2 . 



IV. PROVIDING A REGULAR CENTER AND 
MATCHING WITH A SCHWARZSCHILD 
EXTERIOR 

By looking at Q , it is evident that p = p NFW diverges 
as r — ► 0. A careless examination, from a full general 
relativistic point of view, of the spherical spacetime given 
by JTSJ-lJUJ with p — p NFW , would yield a curvature 
singularity marked by r = 0, associated with the blowing 
up of the Ricci scalar 



R 



8ttG 



pnfw c i 2 r 



p 



(28) 



However, this situation does not apply to NFW halos, 
not only because they are Newtonian systems that must 
be examined within the framework of a Newtonian limit 
of a weak field approach, but because the NFW mass 
density profile is an empirical fitting formula that 
emerges from spherically averaging numerical simulations 
that cannot resolve distances smaller than about 1 % 
of the actual physical radius of the halo |2^|. Hence, 
astrophysicists using this density profile do not actually 
assume infinite central densities, but regard this blowing 
up of p N Fw as an undesired effect due to the extrapolation 
of a fitting formula which (within the resolution limits of 
numerical simulations) provides a rough illustration of 
the fact that density becomes "cuspy" along the central 
halo region, i.e. /? NFW ~ 1/x for 

A practical way to get rid of this inconvenient feature 
is to "replace" a small spherical region < x < xq of the 
NFW spacetime with an "inner" fluid region containing 
the world-line of a regular center. Using the definitions 
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(0)-©, the radius ro of this inner region in terms of M v 
is given by 



= 0.272 x 
kpc V M Q 



0.273 



Xq. 



(29) 



Hence, for halos in the observed range 10 8 M Q < M vir < 
1O 15 M , the choice x = 0.0001 yields 4pc < r < 
340 pc, a very small radius in relation to the virial radii 
of these halos. Thus, since this length scale is much 
smaller than the maximal resolution of numerical sim- 
ulations |23j |. we are able to provide a regular center for 
the NFW spacetime but this does not prevent us from 
studying the effects of its steep density profile in the cen- 
tral region. 

The simplest choice of a spacetime geometry for the 
inner region is a section of a Schwarzschild interior solu- 
tion characterized by the metric (|15|l with 



exp (^j = a o-b v / l _ 



M 



4np c 3 
r 

3c 2 ' 



with Ko = 8 7T G /i c /c and 

/i = p c = const. 



P 



p c 3 6 V 1 - k r' 2 



3 ao — ba Vl — kq r' 2 



(30) 

(31) 
(32) 



where the constants ao, 60 and p c must be selected so 
that this region can be suitably "glued" to the NFW 
spacetime occupying x > xq. 

As we mentioned before, it is customary to disregard 
the asymptotic behavior of NFW profiles, since the virial 
radius is considered to be the physical cut-off radius of 
NFW halos. However, we can construct asymptotically 
well behaved NFW configurations for which p, P, M/r 
and $ tend to zero as x — > 00 (though M most certainly 
will diverge in this limit, since the Newtonian M NFW in 
@ already does). A finite M as r — > 00 can be achieved 
if we match the NFW spacetime at a convenient cut-off 
scale to a Schwarzschild vacuum exterior characterized 
by p = P = and by ljT5|> with 



exp 



2$ 



1 



2GM 



M 



M , (33) 



where Mq is the constant "Schwarzschild mass" . 

Necessary conditions for a smooth matching between 
spacetime regions are given by Darmois matching condi- 
tions [33j , requiring continuity of the induced metric and 
extrinsic curvature of the matching hypersurface 



i a b = n a n b - g ab , 



K„ 



-K Kn c - d , 



(34) 



where n a is a unit vector normal to this hypersurface. 
Since the NFW spacetime must be matched, either to 
(|30|l or to l|33|) . at hypersurfaces marked by constant r, 



we have n a = y/9rr $ r a , hence (|34ll imply that gtt, g' a an d 
g rr (but not necessarily g' rr ) must be continuous at the 
matching hypersurface. Considering (|19[1 - H2HI . this im- 
plies continuity at the matching hypersurface of M, $ 
and P, but not of fi or the anisotropic pressure defined 
in terms of T by iffl)) . 



A. Matching with the inner region. 

It is convenient to assume Ij25(l to be valid at r = 0, 
so that we can characterize the Schwarzschild interior 
solution by 



/i c = p c c 2 + (3/2) P c 



(35) 



Following (|12[) . we can define a characteristic velocity 



K = 4nGp c r 2 s 



so that 



P 



ScPcV c 2 , £ 



Pc C z 



V 2 
6c -§-, 
c 2 



(36) 



(37) 



where 5 C is an arbitrary constant, so that central velocity 
dispersion is a 2 — 5 C V 2 . Hence, for < x < xq we have 



m = fpc 



(38) 



while, for the time being, we assume also T = 0, though 
a nonzero T can be considered for the inner region in the 
case of anisotropic pressure (see section VI-B). Since we 
are considering x < xq ~ 0.0001 <C 1, suitable expres- 
sions for the remaining variables in this region are found 
by expressing ao, 60 i n terms of the parameters in l|35() - 
(|36|l and expanding (|30|l and l|32|) up to first order in x 2 , 
leading to 



P 

V 2 



<*>C + \V 2 (l+lt) <' 



PcV 2 



(39) 

(40) 
(41) 



Following the matching conditions l)34[l. the constants 
p c , V c , $c and S c must be selected so that M(xq), $>(xq) 
and P(xo) continuously match the NFW functions M, $ 
and P at xo . Although, (|3"4"}l do not require this continu- 
ity for p and T, we will still assume it in order to avoid 
an unphysical jump discontinuity of these variables at 
x — xq, as well as all state variables. 



6 



B. Matching with a vacuum exterior. 

A smooth matching with a Schwarzschild exterior at a 
cut-off radius r = n, based on (|34[1 requires 



M(r b ) = M , 
P{n) = 0, 



a 2*(r b )/c 2 



2 G M 
n 



(42) 
(43) 



but do not require /i or T to vanish at rj,. However, a 
jump discontinuity of these variables at an interface with 
a vacuum exterior is much more acceptable than in the in- 
terface between two non-vacuum regions. As we discuss 
in section VI, a convenient cut-off scale for a NFW space- 
time is the virial radius = r vir , so that we can identify 
Mq with M vir . Though, because of the matching with the 
inner region, an r vir selected by means of condition i|43|) 
will not yield (even in the Newtonian limit) Mq — M vir 
with M vir given by JHJ). However, for sufficiently small 
io < 1 the resulting Mq will be approximately equal to 
M vir . 



so that in the limit e — > we recover the Newtonian equa- 
tions 0, © and ifTBf. The system |IS )) -l|I3 f can be inte- 
grated by demanding that M. and V comply with appro- 
priate boundary and initial conditions, so that the NFW 
halo can be smoothly matched with the Schwarzschild 
interior at x = xq and the Schwarzschild exterior at 
r = r vir . Since we have to use the explicit form of Y 
in 144|l . then the analytic or numerical solutions of l|48|> - 
(|49|l for specific choices of T, boundary conditions depend 
on M vir through the definitions and ©. 

The metric function M = V 2 r s M follows from l|4^|> . 
while $ = $ c + V 2 ^> can be obtained by integrating 



d* M + eVx 3 
dx x[x — 2sM]' 



(51) 



The relativistic generalization of the Newtonian rotation 
velocity profile are the velocities of test observers along 
circular geodesies. These velocities are |3,l2i Kot = 
which in terms of the dimcnsionlcss variables becomes 



V 2 = 



v o 



M+eVx 3 
x -2eM 



(52) 



V. POST-NEWTONIAN NFW HALOS 

In order to explore the post-Newtonian limit for the 
system (I2b|) (I27|l , it is useful to work with dimensionless 
variables by rescaling all variables in terms of quantities 
defined at the scale radius r s . 



A. The region x > xo 

Convenient rescalings follow as 
Y 



1 



M 
V 



5 p x[l+x] 2 
M - 3 



Co A M 



ATrdoPorl 3 6 M V 
P 



SoPoV 21 
V 2 ' 



(44) 
(45) 
(46) 
(47) 



with Vo defined in 1|12[1. transforming l|26|l and l|27|) into 



dM 

da; 

dV 

dx 



Y+2 eV ( 48 ) 
[y+feP] [M+eVx 3 ] 2T ^ 



where 



x [x - 2eM] 



v; 



(49) 



(50) 



Since Vo for typical galactic halos ranges from a few 
km/sec to ~ 1500 km/sec, the post-Newtonian cor- 
rections of order V 2 /c 2 will be very small: between 
0(e) ~ 10" 9 and 0(e) ~ 10~ 6 . The post-Newtonian 
system associated with (14811 I49|) can be given as 



dM 

da; 
dV 

da; 



Yx 2 + O(e), 
Y M 2T 



V + 0(e), 



with 



df M 

— = ^2+ 0(e), 



dx 
V' 



= ^ + 0(s) 
x 



(53) 
(54) 

(55) 
(56) 



B. The region < x < xo 

The variables defined in the previous subsection must 
glue continuously at a;o with the interior Schwarzschild 
variables H35(l - (|40|l . Normalizing these variables with the 
same factors as in (I44[) ~(I46 |1 . we have in the region x < 

x < Xq 



Y = Yn = 



M = -Y 



SoPo 

1 + §< 



S c e 



x 3 , 



1 



a; (1 + a; ) ; 



(57) 



V 



Y~ 



(58) 



(59) 
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* « iy (i + f{)^ 



3 u V 2 

From (|57|l and bearing in mind that iq C 1, we have 



(60) 
(61) 



x 



(62) 



implying that £ ~ £ if <5 C ~ xq- Since e is very small 
we can also assume that £ -C 1, so that post-Newtonian 
expressions follow by taking only terms up to O(0 : 



= 




-o(e), 


(63) 


V = 


Y 2 




1 2" 

6 X . 


+ o(0, 


(64) 


* = 








(65) 


V 2 = 




2 z 2 - 


1-0(0- 


(66) 



We examine analytic solutions of the post-Newtonian 
system (I53 fl - H54[) that match continuously with l|63|l - 



VI. ANALYTIC SOLUTIONS. 

For whatever choice of Y and restrictions on P, equa- 
tions and lfS3|l can be integrated, yielding M and $ 
so that A4(xq) and ^(xq) match Q63JI and l|65(l at x = xq. 
Denoting the inner and NFW regions as 



< x < Xq, 
x > Xq, 



(I) , 

(II) . 



(67) 



We have then the following solutions up to orders 0(e) 
0(0: 



M ( i) 
■M(n) 



3x (l + Xq) 2 

(3 + 4x )xq 
3(l + x )2 4 



hi 



o(0= 

1 + x 



1 + Xq 1 + X 



0(e), 
(68) 



logioY 



(a) 



+ 

+ 




+ 




+ 

+ 




X 




(I) 


(II) \ 





1 1 1 1 1 1 1 1 1 

6 -5 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-4 -3 -2-10 1 

lO£inX 


log 




(b) 


0- 






-2- 






-4- 








(I) 


/ (ID 


-6- 






-8- 


+/ 




-icC 

-12l 


+ / 
+ / 
+ / 
+ / 
+ / 




-14- 


/ x 
1 1 1 1 


= x 

logioX 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 



logm 1 ? 

= 



-2 

(C) 









: (i) i y 

+ y 
i+ / 
4-1/ 

+ / 
+ / 
+ / ] 
+ /1 
+ / 1 
- ■+ / 1 


(II) 




- / XtXo 

1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 


logioX 



-2 



(ID 



6a;o(l + Xq) 2 

2 + 3 xq 1 
2(l + x ) 2 + x 
+0(e), 



o(0, 



In 



1 + a (3 + 4i )io 
l + xo 3(l + a;o) 2 



(69) 



Therefore, irrespective of the choice of T and/or assump- 
tions on P, the metric elements for all NFW halo space- 



FIG. 1: Solid curves in panels (a), (b) and (b) respectively 
depict the logarithmic plots of Y, A4 and ^ given by 1441 . 
CT . and in the inner (I) and outer (II) regions 

separated by :ro = 0.0001. In each panel the "pure" NFW 
case without inner region is shown by the curves with crosses. 
Notice how this latter case is practically identical to Y, M 
and ^ in the outer region. 
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times are up to order e 2 



suitable boundary conditions. 



gtt 



s 2* c /c 2 e 2e* = 1 



2eM 



2$ r 



2e* + 0(e 2 ), 



(70) 

l+ 2 -^ + 0(e% (71) 
x 



where the metric functions M. and VP are given by i|68|) 
and l|69[l in the regions (I) and (II) , and we have obtained 
M and $ from M and * by means of |@SJ and (JHOJ. 
Figures la, lb and lc display the normalized density Y 
and the metric functions M. and W in regions (I) and (II) . 

Notice that all functions defined so far reduce to their 
Newtonian NFW forms, as given in section II, in the 
limits xo — > and £ — > 0. Also, while all NFW halos have 
the same rest-mass density Y, the form for the pressure 
depends on the assumptions one might make about T and 



A. Isotropic case 

For r = we have P = P r = P± and so pressure is 
isotropic. In collisionless systems this implies an isotropic 
distribution of velocity dispersion. In this case, l|54|) 
yields the following analytic solution 



V, 



5 C - x 2 /6 

' = X 2 {l+X f 

S c - x\l<o 



+ o(0 

+ T-T + O(e), (72) 



(H) " ^(l + *o) 4 
where Tq is the evaluation at x = xq of the function 



T = 



a 



T(x ,x) = | [ln(l + x)} 2 



[A x - a ] ln(l + x) + [a - 7/2] In a; + 3Li 2 (1 + x) + a Q B 1 + C x , 



(73) 



31n(l + afo) + 



3(l + x )-x 2 
(l + x ) 2 ' 



Ai 



1- 3a; 



lx i 



2x 2 (l + x) 



Bi 



1 - 3a; - 6x 2 
2x 2 (l + x) 



1 - 3a; - 18a; 2 - 13x d 
2a; 2 (1 + x) 2 



and the dilogarithmic function is defined as 

rV hit At 



Li2(y) 



l - 1 



This form of V matches continuously the regions (I) and 
(II). Since the limit of !F{x, xq) as x — > oo depends explic- 
itly on a;o and S c , we need to find appropriate relations 
S c = S c (xq) in order to determine (together with 
H7H ) the asymptotic behavior of the NFW spacetime. 



1. Asymptotically flat configuration 



by 



The asymptotic behavior (a; » 1) of l(72 jl -l|75 )l is given 



V, 



(ii) 



x 2 /6 



X 2 (1 +x )^ 
1 - (3/4) a 



-Fa- 



71" 

T 



4 In x 



16 x 4 



O 



In a; 

x 5 



(74) 



Since in region (II) we have: p — > as x — > oo, an 
asymptotically flat NFW configuration without cut-off 
scales requires that 'P(ii) — > as x — > oo, thus the zero 



order term in (|74|) must vanish, leading to 



l + 6(l + a; ) 4 









■t)] 





7T~ 

T 



17 

t 



in 



,(75) 



where we have used (|73|l and the fact that xq -C 1 in 
order to get this leading term expansion on Xq. Hence, 
(|62[1 implies that £ is smaller than e but of the same order 
of magnitude, so that 0(£) ~ 0(e). The asymptotic limit 
given by lJ73 jl -lfH) )l also implies 9(n) — > and, even if 
M(u) diverges, we have M./x — * h\{x)/x — > so that 
asymptotically we have also g rr — > 1. 



2. Matching with a vacuum exterior. 

For a galactic halo structure the cut-off scale for a 
matching with a Schwarzschild exterior should be the 
virial radius r — r vir (equivalently, x = co), hence condi- 
tion (|4^|) implies 



P(ii)(co) 
^(ii) 



= 0, 

n 

6 

= J" 



f [1 + 6(1 



■ x ) 4 (^(c ) - ^ )] 



(76) 



where T(c ) is F(x , x) evaluated at x = c . For h = 0.7 
and bearing in mind that for all virialized galactic halo 
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logioP 



(a) 















c- 








X 


Co 


-1- 












_2_ 












-3- 


(I) 




/ T T \ 

(II) 






-4- 












-5- 












-6- 


X 


— A() 




S 


\a 












-V- 




1 1 1 1 

6 


1 1 1 

4 


1 1 1 1 1 

-2 


1 1 1 1 






logioP 




(b) 


















o: 








X 


Co 


-i: 
























-3: 


(I) 




(ii) 






-4- 












-5- 


X 


= X 








_6 : 


1 1 1 1 


1 1 1 


1 1 1 1 1 


S 

1 1 1 1 


V 

1 



log 10 x 2 



FIG. 2: Logarithmic plot of normalized pressure V in regions 
(I) and (II) with xo = 0.0001. Panels (a) and (b) respectively 
denote the isotropic (eqs — ) an d anisotropic (eqs l|82[l - 
1841 1 ) cases. The letters A and S depict the asymptotically flat 
case that scales as In x/x 4 and the case with a matching with 
a Schwarzschild exterior at the virial radius (x = Co — 15) so 
that 'P(co) = 0. The dotted curves in the inner regions are 
the NFW cases without an inner region. In the anisotropic 
case we have chosen j3 = 20. The curves deviating from A 
correspond to cases that scale asymptotically as x~ 2 



structures we have 10 8 < M vir /M & < 10 15 , numerical 
values of Cq given by (JSJ) fall in the range 6 



< 



Co 



< 30. 



Thus, considering the function T in l|73(l we can expand 
S c as in J7SJ), leading to 



hi 



0.73 



x In 



(77) 



so that S c is of the same order of magnitude as in H77fl . 
Conditions (|42|1 for the post-Newtonian metric functions 
<EHI>— G3> are given by 



£*(c ) 



M(cq) 



CO 



.M(co) 



(78) 



which imply 



M(r vir ) = A/ v 



2 + 3x 



.2(1 



- a^o) 2 
3£ 
A C 3 



2, 14 
f — In — 

co 1 + 

ln(l + x ) 



Cq _ 1 
x 1 + c 
_ (3 + 4x 2 )x 
x ) 2 

r 



3(1 



+ 0(s 2 ) 



0{e) 



2 1n(l + c ) 1 



co 



1 + co 



(79) 
(80) 



where we have used the definitions of 80 : M viI in Q and ner region, M(r vir ) = M vlr does not hold exactly in the 
©• Notice that, because of the matching with an in- 
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Newtonian limit, though it holds for a very good ap- 
proximation if xq is sufficiently small. We show in fig- 
ure 2a the logarithmic plot of V(x) for the two cases 
considered above (asymptotically flat and matched to a 
Schwarzschild exterior). 



3. Polytropic equation of state. 

The complexity of the expressions in 172(1 and i|73|l do 
not allow us to find out, at first glance, the type of rela- 
tion between P and p. Though, by looking at (|44jl and 
(|74"Jl . the asymptotic behavior V ~ ln(x)/x 4 ~ 1/x 4 and 

Y ~ 1/x 3 indicates a sort of power law relation between 

V and Y that (at least asymptotically) might be simi- 
lar to the polytropic relation 124|) . In order to examine 
the functional relation between Y and V, we provide in 
figure 3a the logarithmic plot of V vs Y (or equivalently 
InP vs pV 2 ), for the asymptotically flat case and the 
case matched with a Schwarzschild exterior, using the 
numerical values xo = 0.1, cq = 8. For theoretical ref- 
erence we show the curve corresponding to a polytropic 
relation l|24|l with n — 10. As shown by the figure, the 
asymptotically flat NFW configuration fits very well this 
polytrope, except for high density values corresponding 
to smaller x. This behavior is reasonable, since closer 
to the center (x close to xo) the NFW density profile be- 
comes cuspy, while polytropic density profiles are charac- 
terized by a "flat core" . In the case of a matching with a 
Schwarzschild exterior, the fitting with a polytrope also 
fails near the boundary x = cq (or r — r vir ), which is 
expected since we have V(cq) = while Y(cq) > 0. 



for Y, M. and \& in regions (I) and (II) that we derived 
in previous sections (ie all equations (|35[) -H71 [) and l|78|) - 
(jEOJ), except for flO")!. (JHHI and ffity ) remain valid for the 
anisotropic case, regardless of the form we might assume 
for r. However, equation (|54|l does involve T and so it 
must be integrated for both regions. 

A useful expression for the anisotropy factor T is the 
ansatz proposed by Ostipkov and Merritt |3J| 



r = 



(81) 



where = rp /r s — cq rp /r vir marks the length scale (nor- 
malized by r vir ) in which the velocities of collisionless 
particles pass from an isotropic regime near x = to a 
radially dominant mode, since T — > 1 (or P± — > 0) as 
x becomes larger. Numerical simulations suggest that 
P±/P — ► 0.6 — 0.8 at about the virial radius x = cq, 
hence we can set ~ k cq with 1.2 < k < 2. 

Although Darmois matching conditions l|34|) allow for 
jump discontinuities of T, we will assume the anisotropy 
factor l|81|) to be continuous at xo and to hold also in 
the domain of the inner region < x < xo with constant 
density. Under this assumption, the form equivalent to 



(i) 



fZS} is 



(3 2 5 c 



x 2 (x 2 + 2/3 2 ) 



x2(l + x ) 4 (^ 2 +/3 2 ) 



+ O(0, (82) 



while in the region (II) the form of V that follows from 
the integration of l|54|l for (|81|l and matches continuously 
with is 



B. An anisotropic example 

Since V is decoupled from M. and \& in the post- 
Newtonian field equations Q53 [l -I|55| ) . all the expressions 



^(ii) = 



(x 2 + f3 2 )V + Q-Q 



x 2 +f3 2 



+ O(e), (83) 



where V Q = V^(x Q ), Q = Q{/3,x ) with 



Q(J3, xo, x) = Y [ln(l + x)] 2 + [A 2 ~ fa 7 o] ln(l + x) + (f3 2 7o - 3 ) ln(x) + /3 2 Li 2 (l + x) - 7o B 2 - C 2 



(84) 



02 

D 2 



1 + 3/3 2 , 3 = 1+2/3 2 , 

2 2 x 2 +0 2 (3x- 1) 



7o = 



3(l + x )-x 2 



+ ln(l + x ), A 2 = 



(70 2 + 2)x 3 + 2 (x 2 - 3x + l) 



2x 2 (l + x) 



C 2 = 



3(l + x ) 2 

(4 + 13/3 2 ) x 3 + (5 + 180 2 ) x 2 + 2 (3x - 1) 
2x 2 (l+x) 2 



2x 2 (1 + x) 



Just as in the isotropic case, we examine the asymptotic configurations, 
behavior of the NFW halos characterized by lIH^ - lp^jl . 
As mentioned before, the forms for M. and \& and the 
metric functions given in i|68 |) -i|7 Hl are valid for these 



11 



1. Asymptotically flat cases. 



As oposcd to Vni) given by (J73J), from ^SJ-JHH wc 



(II) 

have: 'P(ii) — > as x — > oo for any value we might choose 
for the parameters j3, S c and Xo- Hence, all NFW con- 
figurations characterized by the Ostipkov-Merritt ansatz 
(|81|l for r are asymptotically flat. However, by looking 
at the asymptotic behavior of Q 



(ii) 



Co 



9i. 

X 2 

+0 



In x - ln(l + x ) - 2(3 2 C - C x 



2x 4 



In x 

x 5 



(85) 



(^+/3 2 )P -Qo-— p 2 



Ci 



5 x\ 



6(l + x ) 2 



it is evident that the asymptotic behavior depends on Co. 
If Co > 0, then V(u) > decays asymptotically to zero as 
1/x 2 , this case is shown by unlabeled solid curves in fig- 
ures 2b and 3b. However, this case is unphysical because 
(fromJ2s30 the velocity dispersion scales asymptotically 
as a 2 — » Co x and diverges as x — > oo. If we want a 2 — > 
asymptotically, then we must choose Co = 0, leading to 
the same asymptotic scaling "Prm ~ In (a;) /a; 4 as in the 
isotropic case. From (|52l . This case corresponds to the 
choice 



InV 



(a) 



-4- 



-12- 



-1& 



X = iCo 


X 


= x 


/ (ID 




(I) 


X 







InP 



-12- 



-1& 



(b) 



-10 -5 5 10 



In Y 













X = 


CO 




X 


= Xo 






(II) 




(I) 


/ 


* / 


















A 


I s 









-10 



10 lnY 



6c ~ P 2 



■2(3 



12 



(l + x ) 4 



Qo + -r 
6 



, (86) 



and is marked by the letter A in figures 2b and 3b, while 
the curves without mark in these figures correspond to 
various values of Cq > 0. 



2. Matching with a Schwarzschild exterior 

As in the isotropic case, we assume the matching inter- 
face to be r vir so that x = cq. The matching conditions 
(g2J) are given by ffify. leading also to (O and (JHOJ) - 
However, (I43H in the form "P(n)(co) = now implies 



6 = — 

° c ~ I3 2 



2 r^2 



2(3 2 



12 



+ (l + x ) 4 (Qo-Q(c )) 



(87) 



where Q(c ) is given by (|51Jl evaluated at x — c . The 
form of V corresponding to this case is shown as the curve 
is marked by the letter S in figures 2b and 3b. 



3. Polytropic equation of state 

Since "P(n) with Co = follows the same asymptotic 
scaling as in the isotropic case, it is not surprising to find 
that V and Y follow the same approximately polytropic 
relation. However, in the case Co > we see an asymp- 
totic relation of the form P ~ p 2 / 3 usually for x <§C c far 
away from the virial radius (see figure 3b). 



FIG. 3: Plot of In V vs In Y, equivalent to plotting In P 
vs In pV$. Panels (a) and (b) denote the isotropic and 
anisotropic cases for the same parameters as in Fig 2. Let- 
ters S and A identify the curves for the case matched to a 
Schwarzschild exterior and the case that scales asymptotically 
as (In x)/x 4 . The other curves in (b) mark the cases that scale 
as x~ 2 . The dotted curve in the inner region is the "pure" 
NFW without inner region. As a comparison we show a line 
with slope 1.1 (thick grey line) that would correspond to the 
polytropic relation P oc p 1+ll/1 °. Notice how the variables 
characteristic of the NFW profile decaying as (In x)/x 4 ap- 
proximately fit this relation, except near the center and near 
the Schwarzschild matching interface. 



VII. DISCUSSION AND CONCLUSION 

In the previous sections we have constructed ade- 
quate post-Newtonian generalizations for the galactic 
halo models that emerge from the well known NFW nu- 
merical simulations. We have shown how the issues of 
lack of a regular center (because of interpolating an em- 
piric density profile) and an unbounded halo mass can be 
resolved by suitable matchings with a section of an in- 
terior Schwarzschild solution with constant density, and 
with a vacuum Schwarzschild exterior. Even if galac- 
tic halos are essentially Newtonian systems, we feel it 
is important for relativists to see how they can also be 
described and studied in General Relativity within the 
framework of a post-Newtonian weak field regime. Such 
a description can be very valuable in studying their in- 
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teraction with physical effects (gravitational lenses and 
gravitational waves) and dark energy sources, all of which 
lack an adequate Newtonian description. 

Following our proposal that NFW halos satisfy the 
ideal gas type of equation of state (|22|l . we have shown 
cmpiricaly (see figure 3) that outside their central core re- 
gion these halos approximately satisfy the polytropic re- 
lation l|24|l with n « 10. This might be quite significant, 
since virialized self-gravitating systems are characterized 
by non-extensive forms of energy and entropy pH |22| , 
and as mentioned before, stellar polytropes are the equi- 
librium state associated with the non-extensive entr opy 
functional in Tsallis' formalism [MEHHE! (see 
for a critical appraisal). However, the consequences of 
this rough polytropic relation should be looked carefully, 
since stellar polytropes are solutions of Vlassov equation 
with an isotropic velocity distribution 0], while NFW 
halos follow from numerical simulations and exhibit (in 
general) anisotropic velocity distributions (even if these 
anisotropies are not too large Q). In the application 
of Tsallis formalism to self-gravitating collisionless sys- 
tems [13, H3, the free parameter q = (2n— l)/(2n— 3) de- 
notes the departure from the extensive Boltzmann-Gibbs 
entropy associated with the isothermal sphere (which fol- 
lows as the limiting case n — > oo, or equivalently, as 
q — > 1). Assuming Tsallis theory to be correct, the em- 
piric verification provided by figure 3 might indicate that 
in the region outside the central core NFW numerical 
simulations yield self-gravitating configurations that ap- 
proach an equilibrium state characterized by the Tsallis 
parameter q w 1.1. 

While the central cusps in the density profile predicted 
by NFW simulations seem to be at odds with observa- 



tions 0, 0, 0, 0] , there is no conflict between these 
observations and the 1/x 3 scaling of the NFW density 
profile outside the core region. Although the issue of the 
cuspy cores is still controversial, if it turns out that galac- 
tic halos do exhibit flat density cores, their density pro- 
files could adjusted to stellar polytropes and this might 
be helpful in providing a better empirical verification of 
Tsallis' formalism. However, this idea must be handled 
with due case, since stellar polytropes are very idealized 
configurations. 

Although we have only dealt with NFW halos, the 
methodology that we have followed here can be applied, 
in principle, to any Newtonian model of galactic halos. 
For a deeper study of galactic halo models (NFW, as well 
as other empiric or theoretical models), it is important 
to consider a wider theoretical framework, not only using 
a post-Newtonian approach, but including also the usual 
thermodynamics of self-gravitation systems 0, , as 
well as alternative approaches such as Tsallis' formal- 
ism HSilil. This study might provide inter- 
esting theoretical clues for understanding the Statistical 
Mechanics associated with numerical simulations and/or 
gravitational clustering. An improvement and extension 
of the present study of NFW halos are being pursued 
elsewhere |35j j . 
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